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ABSTRACT 

The rich oscillation spectra determined for the two stars, v Eridani and 12 Lacertae, 
present an interesting challenge to stellar modelling. The stars are hybrid objects 
showing a number of modes at frequencies typical for /3 Cep stars but also one mode at 
frequency typical for SPB stars. We construct seismic models of these stars considering 
uncertainties in opacity and element distribution. We also present estimate of the 
interior rotation rate and address the matter of mode excitation. 

We use both the OP and OPAL opacity data and find significant difference in the 
results. Uncertainty in these data remains a major obstacle in precise modelling of the 
objects and, in particular, in estimating the overshooting distance. We find evidence 
for significant rotation rate increase between envelope and core in the two stars. 

Instability of low-frequency g-modes was found in seismic models of v Eri built 
with the OP data, but at frequencies higher than those measured in the star. No 
such instability was found in models of 12 Lac. We do not have yet a satisfactory 
explanation for low frequency modes. Some enhancement of opacity in the driving 
zone is required but we argue that it cannot be achieved by the iron accumulation, as 
it has been proposed. 

Key words: stars: variables: other - stars: early-type - stars: oscillations - stars: in- 
dividual: v Eridani - stars: individual: 12 Lacertae - stars: convection - stars: rotation 
- stars: opacity 



1 INTRODUCTION 

Transport of chemical elements is still not well under- 
stood aspect of stellar interior physics. In the case of the 
upper main sequence there are uncertainties regarding the 
extent of mixing of the nuclear reaction products beyond 
the convective core, known as the overshooting problem, 
as well as, regarding the survival of element stratification 
in outer layers caused by selective radiation pressure and 
diffusion. Closely related is another difficult problem of the 
angular momentum transport because there is a role of rota- 
tion in element mixing. The upper main sequence pulsators, 
P Cephei stars, are potential source of constraints on mod- 
elling the transport processes massive stars. Recently, Miglio 
et al. (2007b) and Montalban et al. (2007) studied sensitivity 
of the oscillation frequencies in B stars to the rotationally 



E-mail: wd@astrouw.edu.pl 



induced mixing. A comprehensive survey of properties of 
these variables was published by Stankov & Handler (2005) . 

Pulsation encountered in f3 Cep stars frequently consists 
in excitation of a number of modes, which differ in their 
probing properties. They are found most often in evolved 
objects, where low order nonradial modes have mixed char- 
acter: acoustic-type in the envelope and gravity-type in 
the deep radiative interior. Frequencies of such modes are 
very sensitive to the extent of overshooting. Their rota- 
tional splitting is a source of information about the deep 
interior rotation rate. Pulsation is found both in very slow 
(< 10 km s _1 ) and very rapid (> 200 km s _1 ) rotating stars. 
Thus, there is a prospect for disentangling the effect of rota- 
tion in element mixing. First limits on convective overshoot- 
ing and some measures of differential rotation have been 
already derived from data on the (3 Cep stars HD 129929 
(Aerts et al. 2003), v Eri (Pamyatnykh et al. 2004 (PHD), 
and Ausselooss et al. 2004). Unfortunately, both objects are 
very slow rotators. To disentangle the role of rotation in el- 
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ement mixing, we need corresponding constraints for more 
rapidly rotating objects. 

It seems that we understand quite well the driving ef- 
fect responsible for mode excitation in /3 Cep stars. The 
effect arises due to the metal (mainly iron) opacity bump 
at temperature near 200 000 K. Yet, the seismic models of 
v Eri, which reproduce exactly the measured frequencies of 
the dominant modes, predict instability in a much narrower 
frequency range than observed. As a possible solution of 
this discrepancy, PHD proposed accumulation of iron in the 
bump zone caused by radiation pressure, following solution 
of the driving problem for sdB pulsators (Charpinet at al. 
1996). PHD have not conducted calculations of the abun- 
dance evolution in the outer layers. Instead, they adopted 
an ad hoc factor 4 iron enhancement showing that it leads to 
a considerable increase of the instability range but also al- 
lows to fit the frequency of the troublesome P2 dipole mode. 
However, recent calculations of the chemical evolution made 
by Bourge et al. (2007) showed that the iron accumulation 
in the j3 Cep proceeds in a different manner than in the sdB 
stars. Unlike in latter objects, in photospheres of /3 Cep stars 
the iron abundance is significantly enhanced for as long as it 
is enhanced in the bump. In fact, the photospheric enhance- 
ment is always greater (P.-O. Bourge private communica- 
tion). However, neither v Eri nor 12 Lac shows abundance 
anomaly in atmosphere except perhaps for the nitrogen ex- 
cess, [N/O]=0.25 ± 0.3, in the former object (Morel et al. 
2006). So there is no clear evidence for element stratifica- 
tion and we are facing two problems: what is the cause of 
element mixing in outer layers of this star, and how to ex- 
plain excitation of high frequency modes {y > 7 cd _1 ) and 
the isolated very low frequency mode at v = 0.43 cd _1 . 

Unlike most of /3 Cep stars, where frequencies of de- 
tected modes are confined to narrow ranges around funda- 
mental or first overtone of radial pulsation, modes in v Eri 
and 12 Lac are found in wide frequency ranges. Both stars 
are hybrid objects with simultaneously excited low-order p- 
and g-modes (typical for /3 Cep stars) as well high-order g- 
modes (typical for SPB stars). The greatest challenge is to 
explain how the latter modes are driven. Our work focuses 
on these two stars, which are best studied but not the only 
hybrid pulsators. Other such objects are 7 Pegasi, t Her, 
HD 13745, HD 19374 (De Cat et al. 2007). Very recently 
Pigulski & Pojmahski (2008) reported a likely discovery of 
five additional objects. 

In the next section, we compare oscillation spectra and 
provide other basic observational data for v Eri and 12 Lac. 
Seismic models of these two objects are presented in Section 
3, after a brief description our new treatment of the over- 
shooting. In the same section, we give for both objects the 
estimate of rotation rate gradient in the abundance vary- 
ing zone outside core. Section 5 is devoted to the problem 
of mode excitation. Finally, in Section 6, after summarizing 
the results, we discuss the matter of element mixing in outer 
layers and what is needed for progress in B star seismology. 

2 OSCILLATION SPECTRA AND OTHER 
OBSERVATIONAL DATA 

The multisite photometric (Handler et al. 2004, Jerzykiewicz 
et al. 2005) and spectroscopic (Aerts et al. 2004) campaigns 
resulted in detection of nine new modes in the v Eri oscilla- 



tion spectrum, in addition to four, which have been known 
for years. Furthermore, data from these campaigns allowed 
to determine (or at least constrain) the angular degrees for 
most of the modes detected in this star. A schematic oscil- 
lation spectrum for this star is shown in the upper panel 
of Fig. 1. The angular degrees, £, are adopted after De Rid- 
der et al. (2004). In the lower panel of Fig. 1, we show the 
corresponding spectrum for 12 Lac. All data, including the 
^-degrees, are from Handler et al. (2006). The two spectra 
are strikingly similar. However, owing to the appearance of 
two complete I = 1 triplets, the one for v Eri is more reveal- 
ing. 

The mean parameters of the two objects are also sim- 
ilar. In Fig. 2 we show their positions in theoretical H-R 
diagram with the errors based on the parallax and pho- 
tometry data (see caption for details). In the same figure, 
we show also the evolutionary tracks for selected seismic 
models of the stars which we will discuss in the next sec- 
tion. Both objects have nearly standard Population I atmo- 
spheric element abundances. In particular, Niemczura and 
Daszyhska-Daszkiewicz (2005) quote the following values of 
the metal abundance parameter [M/H]: 0.05 ±0.09 for v Eri 
and — 0.20±0.10 for 12 Lac. These numbers agree within the 
errors with earlier determination by Gies & Lambert (1992). 
These authors also provide the T e g values which are left of 
the error box shown in Fig. 2. The adopted by us position in 
the H-R diagram indicate that both stars are in advanced 
core hydrogen burning phase. 

The projected equatorial velocity in both stars is low. 
Abt et al. (2002) give 20 and 30 km s" 1 , for v Eri and 
12 Lac, respectively, with the 9 km s _1 formal error. The 
corresponding numbers quoted by Gies & Lambert (1992) 
are (31 ± 3) km s" 1 and (39 ± 14) km s" 1 . For v Eri, the 
two values barely agree (within the errors) but are in con- 
flict with the seismic determination of PHD, who derive the 
value of about 6 km s _1 for the equatorial velocity which 
we believe to be a more reliable value. In Section 4, we will 
present our refined seismic estimate of the rotation rate for 
this star and the first such determination for 12 Lac. 



3 SEISMIC MODELS AND INTERIOR 
ROTATION RATE 

For specified input physics, element mixing recipe, and 
heavy element mixture, we need the four parameters: mass, 
the age, and the initial abundance parameters, (X,Z), to 
determine stellar model, if dynamical effects of rotation are 
negligible. The latter is very well justified for the two stars 
we are considering in this paper. This might suggest that 
for testing the mixing recipe, frequencies for more than four 
pulsation modes are needed. This is not true because there 
are additional (non-seismic) observational constraints on the 
models. On the other hand, there are other uncertainties of 
modelling. The most important concerns opacity. 

In the approach adopted in this paper, like in PHD, we 
fixed X at 0.7, as the changes within reasonable ranges do 
not have any significant effect on frequencies. Furthermore, 
we considered models in the expansion phase, so that in- 
stead of the age we could use T c g as the model parameter, 



The two hybrid B-type pulsators v Eri and 12 Lac 3 



n 1 r 



n 1 r 



n 1 r 



A, 



10 



10 



i 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



1 



J LI ! I I I I I I — 



-5.6 



5.8 



6.2 



J I L_ 



J I L 



J I — J_ 



uKri 



i i i i i 



i i r 



t i r 



n i r 



~i 1 r 



100 



A 



10 - 



10 



n i i I i i i i I i i i i I i i i r 
1 



2 



;; 2 



J_1_J I I I I I I LU 



1 or 2 



J I I I Ll L 



4.5 



5.5 



J I L 



12 Lac 



1 or 2 



J I I Li 







4 

v, [c/d] 



6 



5 



Figure 1. Oscillation spectra of v Eri and 12 Lac based on Jerzykiewicz et al. (2005) and Handler ct al. (2006) data, respectively. The 
numbers above the bars are the most likely £- values, as inferred from data on amplitudes in four passbands of Stromgrcn photometry. 



which is more convenient. There are two modifications in our 
modelling. Firstly, in order to asses uncertainty in opacity, 
we constructed models employing both OP (Seaton 2005) 
and OPAL (Iglesias & Rogers 1996) data and two heavy 
element mixtures, the traditional one of Grevesse & Noels 
(1993, hereafter GN93) and the new one of Asplund et al. 
(2004, 2005, hereafter A04). Secondly, we allow for partial 
mixing in the overshoot layer. In our reference model, we 
adopted the OP opacities calculated for the A04 mixture 
and assumed no overshooting. 

To describe partial mixing, we use an additional ad- 
justable parameter. In the partially mixed layer, extending 
from the edge of the convective core at the fractional mass q c 



to distance d ov = a ov H p (q c ), the fractional hydrogen abun- 
dance was assumed in the form 

X = X c + (q-q c ) w [a + b{q-q c )]. (1) 

The input parameters are w and a ov ■ The a and b coef- 
ficients in the expression above were calculated from the two 
input parameters by imposing continuity of X and its deriva- 
tive at the top of the layer. The limit w — > oo corresponds 
to the standard treatment. We wanted to see, whether this 
additional degree of freedom may help us to fit frequency 
of a certain troublesome mode. In more general sense, it is 
important to know whether asteroseismology may teach us 
about overshooting more than only its distance. 
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Figure 2. v Eri and 12 Lac in the theoretical H-R diagram. 
The the lcr error box for the first star is similar to that adopted 
in PHD, except that it is moved downwards by 0.1 because we 
do not include now the Lutz-Kelker correction to measured par- 
allax. For the second star, the error box is adopted after Han- 
dler et al. (2006). The evolutionary tracks correspond to seismic 
models (dots) calculated assuming no overshooting and the ini- 
tial hydrogen abundance X = 0.7. The v Eri model has mass 
M = 9.625M0 and the heavy element mass fraction Z = 0.0185, 
and that of 12 Lac has M = 11.55M© and Z = 0.015. 



In our seismic estimate of the interior rotation rate, we 
make a similar assumption as PHD, that is in chemically ho- 
mogeneous envelope we assume a uniform rotation rate, Q e , 
and allow its sharp inward rise in the /i-gradient zone. Here, 
we additionally assume uniform rotation rate, Q c , within 
the convective core and a linear dependence of the rate f2 
on the fractional radius, x, in the adjacent p-gradient zone 
between x = x c and x = x Ct o, which is the fractional radius 
corresponding to convective core mass at ZAMS. Thus, wc 
use the following simple expression for interior rotation rate 



fie 



if x < x c 

if x c < x < x Ct o 

if X > Xc.o 



(2) 



and Q e from measured or in- 



Our aim is to determine fi c 
ferred rotational splitting. 

3.1 v Eridani 

Seismic models of this star constructed by PHD made use 
of the three frequencies: 5.637, 5.763, and 6.224 cd _1 , which 
were associated with (£ = 1, gi), (£ = 0, pi), and (£ — 1, pi) 
modes, respectively. Only the mean frequencies of the i = 1 
triplets are matter for seismic models. Individual frequencies 
probe internal rotation. The frequency of 7.898 cd _1 associ- 
ated with I — 1 degree could not be reproduced with stan- 
dard models. PHD noted that the localized iron enhance- 
ment may remove the discrepancy but, as we pointed out 
in the introduction, this is unlikely. Ausseloos et al. (2004) 
succeeded in fitting the four frequencies only for rather unre- 
alistic chemical composition parameters and effective over- 
shooting. 

In Fig. 3, we show results of exploration of the two- 
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Figure 3. Differences between observed and calculated frequen- 
cies of v Eri. All models have been calculated with the OP opacity 
data and the A04 heavy clement mixture. 



parameter freedom in description of the overshooting for fit- 
ting six measured frequencies. All models are constructed to 
fit the three lowest frequencies. Increase in a ov and/or in w 
means more efficient overshooting and reduction of the O-C 
values for the higher frequency mode. The price to pay for 
the improved fit is shown in Fig. 4. All models calculated 
with Q ov > are outside the lcr error box in the H-R dia- 
gram. Still, we regard efficient overshooting and cooler star 
as a possible way of fitting the three high frequency modes. 
There is, however, another problem regarding this part of 
the spectrum. The two peaks at the highest frequencies can- 
not be interpret in terms of remaining components of the 
I = 1, p2 triplet because the frequency separation should be 
similar to that in the pi case and it is much smaller. We 
checked that there is no other low degree mode in this vicin- 
ity, which could be associated with any of the two peaks. 
Any one of them may be a component of the triplet but 
certainly not both. 

It is important to notice in Fig. 3 the difference between 
modes in the O-C changes. For the two highest frequency 
modes, the shift in w from 2 to 8 has a similar effect to the 
shift in a ov from 0.1 to 0.2. For the mode at 6.732 cd _1 
the effect of the latter shift is much larger and results in 
sign change of O-C. We identified this mode as £ = 4, gi. 
This mode was detected only with the data from the second 
photometric campaign (Jerzykiewicz et al. 2005) and no £- 
value from the amplitude data has been assigned to it. In 
our seismic models there is no lower degree mode near this 
frequency. We may conclude that there is a prospect for 
getting additional constraints on overshooting once we have 
enough frequency data. 

As already noted by PHD, the peak at 7.2 cd~ x may 
only correspond to £ = 2, pi mode. Its O-C values are ob- 
tained assuming m = and they would be much closer to 
zero if m = —1 (retrograde) identification is adopted. At this 
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Figure 4. v Eri in the theoretical H-R diagram. The the 1<t er- 
ror box is inferred from the observational data. The evolutionary 
tracks correspond to seismic models (dots) calculated for different 
values of the overshooting parameters. 

point, however, we are reluctant to use this mode in seismic 
model construction. Higher frequency modes are more sensi- 
tive to uncertainties in outer layers and we know that some 
modification of the structure of this layer is required because 
we have to solve the excitation problem. 

In Fig. 5, we show the effect of choice of opacity data 
and the heavy element mixture on seismic models. We may 
see that the latter is far less important. Models calculated 
with the OPAL data are hotter and therefore may accommo- 
date overshooting and stay within the error box. However, 
relaying on the OPAL data will not help the I = 1 , p2 mode 
frequency fit, as the entries in Table 1 show. The smaller 
difference between measured and model effective tempera- 
ture, AlogT c ff, is always associated with greater O-C for 
this mode. Therefore, it seems that if future measurements 
confirm the current value of the effective temperature, the 
only way of reconciling the seismic and non-seismic observ- 
able is an increase of opacity in outer layers. In PHD this 
was accomplished by an artificial enhancement of the iron 
abundance which had virtually no effect on frequencies of 
the three modes used in the seismic model but caused the 
desired -0.1 cd _1 frequency shift of the I = 1, p2 mode. 

For a specified stellar model, the two parameters, Q c 
and r2 e , describing internal rotation may be directly deter- 
mined from measured frequency splitting, S = 0.5(V+ — ^-), 
of the two £ — 1 modes. To this aim we use the relation 

S= [ dxK.(x)tt/2ir, (3) 
Jo 

where K. is the familiar rotational splitting kernel and fi is 
given by Eq. (2). The kernels for the these modes, which were 
plotted by PHD, differ significantly and this is essential for 
an accurate determination of Q c and Q e . 

There are two differences relative to that work. Firstly, 
we now have measured frequency of the m = — 1 component 
of the higher frequency triplet (pi) and fi c denotes now the 
rate in the convective core and not the mean value over 
the [0, x Ct o] range. In Table 1, we give the the values of the 
f2 c /f2 e ratio on the surface equatorial velocity, Vrot = QeR, 
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Figure 5. Influence of the choices of the opacity data and the 
heavy element mixture on seismic models of v Eri. The evolution- 
ary tracks were calculated assuming a ov = 0. 



Table 1. The O-C (in cd 1 ) for the £=1, P2 mode frequencies and 
the differences between observed effective temperature and that 
of the seismic model, AlogT cff = logT off _ obs - logT offjCal . The 
observational uncertainty in effective temperature is A log T e ff = 
0.011. Models with a ov > were calculated with w = 8. V r ot is 
given in km s — 1 . 



k mixture a OVCT O-C AlogT c ff Q c /Q e V r ot 

OP A04 0.0 -0.127 -0.0103 5.55 5.93 

OP GN93 0.0 -0.151 -0.0075 5.36 5.95 

OPAL A04 0.0 -0.188 -0.0044 5.36 5.99 

OP A04 0.1 -0.085 -0.0159 5.82 5.91 

OP A04 0.2 -0.034 -0.0244 5.78 5.93 



calculated for the selected models. The values differ only lit- 
tle. Certainly, the main uncertainty of our inference is the 
^-dependence of its rate, which was assumed in a very sim- 
plistic form (Eq. (2)). 



3.2 12 Lacertae 

In Table 2 we list frequencies and angular degrees of 
the five dominant modes in the oscillation spectrum of this 
star. The frequencies and the £ values are from Handler 
et al.(2006) analysis of their photometry data. The m val- 
ues were kindly provided to us by Maarten Desmet and are 
based on his analysis of line profile variations. He also pro- 
vided us the value of V e = 52 ± 5 km s _1 for equatorial 
rotation velocity at the surface. 

The first attempt to construct seismic model of 
12 Lac (Dziembowski & Jerzykiewicz 1999) concentrated 
on interpretation of the very nearly equidistant triplet 

/3) fs)- T ne distances to the side peaks are -0.1558 and 
+0.1553 cd -1 . Determination of the £ values by Handler et 
al. (2006) made all proposed interpretations invalid. This is 
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Table 2. The frequencies and mode identifications for 12 Lac. 



I dent . 


Frequency 
(cd- 1 ) 




e 




III 


h 


5.179034 




i 




1 


h 


5.066346 




i 







h 


5.490167 




2 




2 or 1 


h 


5.334357 












h 


4.24062 




2 




? 


fe 


7.40705 


1 


or 


2 


? 


h 


5.30912 


2 


or 


1 


? 


h 


5.2162 


4 


or 


2 


? 


h 


6.7023 




1 




? 


fio 


5.8341 


1 


or 


2 


? 



not a rotationally split triplet. Equidistancy enforced by a 
nonlinear (cubic) resonant mode coupling is also excluded 
by the I values. Thus, just a mere coincidence seems the 
only explanation. 

Construction of seismic model of this star is more com- 
plicated than that of v Eri because it cannot be done sep- 
arately of estimate of the internal rotation rate. Therefore, 
in this paper we limit ourselves to providing evidence for 
a significant inward rise of the rotation rate. In this pre- 
liminary analysis, we consider only models calculated with 
the A04 mixture with Z — 0.015 and without overshooting. 
Upon ignoring all effects of rotation, we determined stellar 
mass and effective temperature by fitting frequencies of the 
I = 0, pi (fundamental) and £ = 1, gi modes to f± and fi, 
respectively. We checked that with the constraints on L and 
T e ff there is no alternative identification of radial orders of 
the two modes and that models in the contraction phase are 
excluded. 

The rotation rate of 12 Lac, as determined from spec- 
troscopy and implied by the £ = 1 mode splitting, is fast 
enough so that the second order effects must be considered 
even at this preliminary phase. Within the O 2 accuracy we 
may write 

u m = u(0) + D () (n 2 ) + Sm + D 2 (n 2 )m 2 , (4) 

where v(0) is the frequency calculated ignoring all effects of 
rotation, S is given in Eq. (3), the Do and D2 terms include 
second order effect of the Coriolis force and the lowest order 
effect of the centrifugal force. In our treatment we rely on the 
formalism described by Dziembowski & Goode (1992) but 
simplified to the case of the shellular (fi = fi(r)) rotation. In 
our treatment the surface-averaged effect of the centrifugal 
force is incorporated in evolutionary models and effects of 
the Coriolis force are included in zeroth-order equation for 
oscillations. So the only rotational frequency shift for nonra- 
dial modes calculated by means of the perturbation theory 
is that caused by centrifugal distortion. For radial modes we 
add only the shift 4/3u 2 ot /uo caused by feedback effects of 
the toroidal displacement induced by the Coriolis force. 

We began assuming an uniform rotation and deter- 
mined its rate from the fi — fi frequency separation using 
Eq. (4) with coefficient calculated in the model reproducing 
the the frequencies of /2 and f± identified as corresponding 
to (n,£,m) = ( — 1,1,0) and (1,0,0) modes. In this model 
the frequency of the (n, £, m) = ( — 1, 2, 1) mode is very close 
to f 3 . However, the inferred equatorial velocity of over 95 
km s" 1 is much larger than that determined from line pro- 
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Figure 6. The splitting kernels for three noradial modes: i = 
l,gl (solid line), £ = 2,gi (dotted line) and £ = 2,g2 (dashed 
line). Triangles mark boundaries of the /i-gradicnt zone. 



Table 3. Models fitting the four dominant mode frequencies (/1 — 
fn) at Sl c = Q e and Q c = A.7Q e . The quantity A5 is the frequency 
mismatch between /s and V-2,2,2 in cd -1 . 



Clc/fle Vrot M/Mq logT cff log L / Lq A 5 TO/ 3 

1.0 95.5 11.771 4.3878 4.194 0.113 1 
4.65 47.1 11.777 4.3891 4.198 -0.134 2 



file changes. An additional problem with this solution is the 
frequency distance of nearly 0.3 cd -1 between /s and the 
nearest quadrupole mode, (n,£,m) — (—2,2,2). We believe 
that these two discrepancies exclude the hypothesis of uni- 
form rotation. 

The splitting kernels for the three nonradial modes plot- 
ted in Fig. 6 are very different. This means that there is a 
potential for deriving significant constraints on differential 
rotation. However, we have direct measurement of only one 
rotational splitting, the others must be inferred through the 
frequency fitting. At the equatorial velocity of rotation of 
about 100 km s _1 , the mean effect of centrifugal force in- 
duces frequency shifts exceeding 0.01 cd" 1 . Thus, if we aim 
at such a precision of the fit, we have to take into account dif- 
ferential rotation in evolutionary models. This we leave for 
our future work. In the present work, the effect of the dif- 
ferential rotation described by Eq. (2) was rigorously treated 
only for calculation of the rotational splitting while the mean 
shifts of the multiplets induced the mean centrifugal force 
were determined assuming uniform rotation corresponding 
to f2 e . 

We considered different values of the Q c /Q e ratio start- 
ing with 1 and used the /1 — /2 frequency difference to deter- 
mine O e . Then, we adjusted the mass and effective temper- 
ature to fit the frequencies of the four dominant modes. As 
we may see in Table 2, there are two possible in values for 
the /3 mode. With m = 1 the fit is achieved at Q c /£l e — 1 
and with m = 2 at £l c /£l e = 0.465. Data on the two seismic 
models are summarized in Table 3. 

We believe that our second model is more realistic be- 
cause its surface equatorial velocity agrees within la with 
observations. The mismatch between /s and calculated fre- 
quency of the (n, £, m) — ( — 2, 2, 2) mode is somewhat higher 
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Figure 7. Frequencies of high frequency modes detected in 12 
Lac (veritical lines) and frequencies off all I ^ 2 modes (dots) in 
the two models (see Table 3). Identified modes arc shown with 
encircled dots. Note than for the model with uniform rotation 
(upper panel), we have no identification for the /8 peak. 

than in the first model. However, it should be stressed that 
this mode is most sensitive to detailed treatment of rotation 
in the inner layers. We may see in Fig. 6 that the splitting 
kernel, which is similar to the mode energy distribution, for 
this mode is very much concentrated in the /i-gradient zone. 
Frequencies of such modes are very sensitive also to over- 
shooting. 

Furthermore, as Fig. 7 shows, the second model repro- 
duces closer frequencies of the remaining modes detected 
in 12 Lac. In this model, fg may be identified as the 
pi, £ = 1, m = — 1 mode, while ft and fg may be, respec- 
tively, the m = and m = — 1 components of the gi, I = 2 
mode. It should be noted that this identification for ft im- 
plies that this mode may be coupled with radial mode asso- 
ciated with f±, which is very close. The slight frequency shift 
caused by this proximity (see e.g. Daszyhska-Daszkicwicz et 
al. 2002) should be included in the fine frequency fitting. 
There is a plausible identification of /io as a pi, I = 1. The 
only troublesome peak is fg but only if indeed £ = 1, as 
Handler ct al.(2006) proposed. There would be no problem 
if it was 1 = 2 because then the peak could correspond to a 
fundamental (n = 0) mode. 



4 MODE INSTABILITY 

Satisfactory seismic models should account not only for mea- 
sured mode frequencies but also for their excitation. In prac- 
tice, it means that modes we want to associate with observed 



frequencies should be unstable. As a measure of mode insta- 
bility we use the normalized growth rate 

- W 
r] - r R \dw\ dr ' 

JO I dr | 

where W is the global work integral. It follows from this 
definition that rj varies in the range (— 1, +1) and it is > 
for unstable modes. The value of r\ is more revealing than 
the usual growth rate, because it is a direct measure of ro- 
bustness of our conclusion regarding mode stability. 

In Fig. 8, we show rj for i = 1 and 2 modes calculated 
in wide frequency ranges for our best models of the two 
stars. The overall pattern of the r)(y) dependence in two 
stars is similar. There is the instability range accounting for 
excitation of most of the observed high frequency modes. 
There is also a maximum of r\ in the low frequency range not 
far from the position of the low frequency modes detected in 
the two objects. There is enough agreement to feel that we 
are on the right track toward satisfactory seismic models. 
However, there is a need for improvement. 

Comparing our results for v Eri with those presented by 
PHD, we find a significant increase of the instability range. 
Now, only the highest frequency modes near v = 7.9 cd _1 
and the mode at v = 0.43 cd -1 fall into stable ranges. This 
change is a consequence of the usage of the OP data. That 
these data lead to wider instability ranges than the OPAL 
data, has been already noted before (see e.g. Miglio et al. 
2007a). With the OP data we have unstable 1 = 2 high-order 
g- modes (n=-14 to -20), but their frequencies are still too 
high by factor at least 1.4. The measured frequency nearly 
coincides with r\ maximum for 1=1 but the value is less than 
zero, ?7 m ax = —0.14. With the OPAL data we find maximum 
value at nearly the same frequency but 7/ m ax = —0.33. Per- 
haps, there is still a room for improvements in opacity that 
would render modes at the lowest frequency and the highest 
frequencies unstable. 

All modes detected in 12 Lac, except of that at v = 
0.355 cd _1 , occur in the unstable frequency range. However, 
the difficulty in explaining excitation of the low frequency 
mode is even higher than in the case of v Eri. In part this is 
due to lower mode frequency and in part to the lower heavy 
element mass fraction adopted [Z = 0.015). 



5 DISCUSSION 

We believe that our seismic models of v Eridani yield a 
good approximation to its internal structure and rotation 
but there is room for improvement and a need for a full ex- 
planation of mode driving. We did not fully succeed in the 
interpretation of the oscillation spectrum of v Eridani with 
our standard evolutionary models and our linear nonadia- 
batic treatment of stellar oscillations. The frequency misfit 
between the high frequency peak to the I = l,p2 mode is 
much reduced in models built with the OP opacity data al- 
lowing large overshooting but such models are much cooler 
(2cr ) than the mean colour of the star implies. Such models 
also nearly avoid the driving problem. Similar conclusions 
regarding models allowing large overshooting distance were 
reached by Ausseloos et al. (2004) . As for the consequences 
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Figure 8. Normalized growth rate as a function of mode fre- 
quency in models of v Eri (top) and 12 Lac (bottom) (the same as 
used in Fig. 2). The frequencies of the detected modes are marked 
with vertical lines starting at the r\ = axis. The length is pro- 
portional to the mode amplitude. 



of usage of the OP instead of OPAL opacity data for mode 
instability, our finding agrees with more general observation 
of Miglio et al. (2007a) that models using the former data 
predict wider frequency ranges. 

An assessment of the inward rise of the angular rate in 
the /i-gradient zone was made for both, v Eri and 12 Lac, 
yielding the values around five for the ratio of the core to 
envelope rate. The surface equatorial velocities in the two 
stars are very different. For the former object our seismic 
estimate gives about 6 km s _1 . For 12 Lac, our estimate of 
about 50 km s _1 is less certain because data on rotational 
splitting are much poorer but the value agrees with estimate 
of Desmet (private communication), based on his analysis of 
line profile changes. There is a large difference in rotation 
rate between the two stars but unfortunately the accuracy 
of our modelling is insufficient for addressing the question of 
the relation between the overshooting distance and rotation. 

The surface rotation rate in v Eri is indeed very low 
and this was the reason why PHD suggested that chemical 
element stratification in this star may be sustained. Hovever, 
a closer look at the problem reveals that the effect of rotation 
on element distribution should not be ignored, even at the 
equatorial velocity of few km s _1 . As Bourge ct al. (2007) 
showed for their 10M Q model, the two-fold excess of the 
iron abundance in the driving zone is produced in the time 
scale comparable with main sequence life time. The process 
is slower than in less massive objects (Seaton 1999) but fast 
enough so that if it goes unimpeded a significant excess of 
iron could be created. The effect of rotation could be ignored 
if the time, r mc , for meridional circulation to travel from the 
depth where most of iron is moved up (T w 4 x 10 K) to the 
bottom of the convective zone (T « 2 x 10 5 K) around the 
iron opacity bump is longer than the star age. To estimate 



r mc we use the well-known expression (see e.g. Tassoul 2000) 
for the speed of the meridional flow 

eR 



V mc = — f(x)P 2 (cos8), 

TKH 



(5) 



where e = Q. 2 R 3 /GM, t kh = g g£- is the Kelvin-Helmholtz 
time, and 



f{x) ~ 3 (ainrX i 



V/Vad 



n 2 \ 

2-wGp ) 



(6) 



The last term, known as the Gratton-Opik term, is kept 
though it is of higher order because it is important near the 
surface where density p is low, even if e << 1 and otherwise 
linearization is valid. In our model of v Eri, the the Gratton- 
Opik term is -3.4 at the depth where most of iron is supposed 
to be pushed up. 

_ tkh Ax 

£ l/l 

We evaluate r mc at the layer where T = 4 x 10°K. and 
for Ax use the distance to the bottom of the convective 
zone. For the v Eri model, we have e = 1.2 x 10~ 4 , t K h = 
5.9 x 10 4 yr, x b = 0.907, Ax = 0.05, / = -15. With this 
numbers, we get from Eq. (6) r mc = 1.14 x 10 6 y, which is 
much less than the age of the star, 1.75 x 10 7 y according to 
our seismic model. Thus, even in this slowly rotating star 
the effect of rotation cannot be ignored. For our best model 
of 12 Lac, we find r mc by about four orders less than in v Eri. 
Yet, as we may see in Fig. 8, the problem with driving the 
low frequency mode is even greater than in the case of v Eri. 

Meridional flow and/or turbulence developed through 
instability induced by rotation are expected to prevent ac- 
cumulation of iron in the driving zone and in photosphere. 
This is a likely explanation of apparently normal chemical 
atmospheric composition of /3 Cephei stars, including such 
slow rotators as v Eri. This is also the reason to reject iron 
accumulation in the driving zone as the solution of the driv- 
ing problem, whose solution must be searched in a different 
way. 

We attach greatest hope for solution of the driving prob- 
lem with improvements in stellar opacity data. Thus, we 
would like to encourage further effort in this field. Reliable 
opacity data are essential for B star seismology. Data on 
mode frequencies in v Eri and 12 Lac are abundant and ac- 
curate enough for probing rotation and structure of the //- 
gradient zone above the convective core. However, to answer 
open questions regarding macroscopic transport of elements 
and angular momentum, we need accurate microscopic input 
data, especially on opacity. Only if we have measurements 
of rotational splitting, our inference on differential rotation 
does not rest on very precise model of stellar structure, but 
such data are rare. Regarding observational work, we view as 
most important new spectroscopic observations of the two 
stars aimed at improving the value of the mean effective 
temperature and leading to the I and m identifications for 
a greater number of modes. 
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